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Abstract 

We review the physics of polarization observables in high energy reactions in 
general and discuss the status of NLO one-loop corrections to these observables 
in specific. Many high order radiative corrections exists for rates but not many 
NLO radiative corrections exist for polarization observables. The radiative correc- 
tion calculations for polarization observables are somewhat more complicated than 
those for rates. They tend to be smaller than those for the rates. In most of the 
examples we discuss we include mass effects which significantly complicate the ra- 
diative correction calculations. We elaborate a general scheme which allows one to 
enumerate the number of independent density matrix elements in a reaction and 
provide explicit examples of angular decay distributions in self-analyzing decays 
that allow one to experimentally measure the density matrix elements. We provide 
examples of reactions where certain density matrix elements are only populated at 
NLO or by mass effects. In our discussion we concentrate on semileptonic bottom 
and top quark decays which are linked to leptonic \i and r decays through a Fierz 
transformation. 
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1 Introduction 



Most of the NLO radiative corrections to rates have been done quite some years ago. In 
fact, two-loop NNLO calculations for rates are now becoming quite standard. In e + e~- 
annihilation one is even pushing hard to determine the NNNNLO corrections to the R 
ratio [1, 2]. Contrary to this, many NLO radiative correction calculations to polarization 
observables involving also massive quarks have only been done in the last few years. 

It is clear that higher order radiative corrections to unpolarized observables will always 
be at the center of attention because data on unpolarized observables are as a rule much 
more accurate than data for polarized observables. Nevertheless, as experimental data on 
polarization observables has been accumulating over the years there is an evolving need 
for radiative corrections to polarization observables. One of the reasons that polarized 
radiative correction were lagging behind unpolarized radiative corrections is that the 
computational effort in the calculation of radiative corrections to polarization observables 
is larger than that for rates. For once, one cannot sum over the spins of intermediate 
states whose polarization one wants to calculate. One therefore cannot make use of the 
powerful unitarity method to calculate radiative corrections from the absorptive parts of 
higher order loop graphs. Further, the definition of polarization observables brings in 
extra momentum factors in the integrands of the requisite phase space integrals which 
makes life more difficult. This is particularly true if the masses of particles in the process 
cannot be neglected. For example, in the process t — > b + W + the longitudinal component 
of the polarization vector of the top quark along the VF-direction is given by 



where the denominator factor \q\ = y qfi — mfy comes in for normalization reasons (q is the 
momentum of the W + ). It is quite clear that such square root factors lead to nontrivial 
complications in the phase space integrations. Similar square root factors appear when 
projecting onto the polarization states of the W + . 

Inclusion of mass effects as e.g. in the semileptonic decay b — > c+l~+vi {m c /m h ps 0.30) 
or in the leptonic decay t~ — > /i~ + v T + (m^/m T ps 0.07) render the analytical 
calculations considerably more complicated. The villain is the Kallen function (m\ + 
m\ + m\ — 2{m\m^ + m\m\ + m^wij)) 1 / 2 which is brought in by extra three-momentum 
factors. In the case that e.g. 1713 — > the Kallen function simplifies to (mi — 1712) 2 which 
leads to an enormous simplification in the phase space integrals. 

Physically speaking mass effects become large in regions of phase space where the 
massive particles become nonrelativistic. For example, in the leptonic decay of the muon 
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fjT — > e _ + z/ M + z/ e , the mass of the electron cannot be neglected in the threshold region 
where the energy of the electron is small. This is illustrated in Fig. 1 where the longitudinal 
polarization of the electron is plotted against the scaled energy x = 2E e /m lx of the electron 
[3]. We have chosen a logarithmic scale for x in order to enhance the threshold region. 
In the threshold region the longitudinal polarization deviates considerably from the naive 
value P l e = — 1. The radiative corrections in the threshold region can be seen to be quite 
large. 
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Figure 1: Longitudinal polarization of the electron in leptonic muon decays at LO and 
NLO as a function of the scaled energy x = 2E e /m fl [3]. The LO curve is very well 
described by the functional behaviour P\ = —(5. 

All our results are given in closed analytic form. In the days of fast numerical comput- 
ers one may rightfully ask what the advantage of having closed form expressions is. One 
of the advantages is that all requisite mass and momentum limits can be taken in analytic 
form. Thus the analytic results can be subjected to tests against known limiting cases 
providing for the necessary checks of the full results. A second issue is user-friendliness. 
The analytic expressions are simple enough to be incorporated in numerical programs by 
the prospective user. Parameter values such as masses and coupling constants can be var- 
ied at will by the user in his own program without having to refer to numerical programs 
written by others. We have checked that our analytical formulas are numerically stable 
even in the small and large mass limits. 

NLO corrections to rates can be quite large. NLO corrections to polarization observ- 
ables have a tendency to be somewhat smaller. The reason is that polarization observables 
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are normalized quantities. They are normalized with regard to the total rate. The numer- 
ators and denominators in the relevant polarization expressions tend to go in the same 
direction. The reason is that the NLO numerators and denominators are dominated by 
the soft gluon or soft photon contributions which are universal in the sense that they mul- 
tiply the relevant Born term expressions and thus cancel out in the ratio. For example, 
for t -> b + W + the 0(a s ) correction to the rate is -8.5% and (-2.5% -=- +3%) for the 
polarization observables [4, 5, 6]. 

At NLO some density matrix elements become populated which vanish at leading 
order. An example is again the decay t — > b + W + where the rate into the right-handed 
W + vanishes (T(W R ) = 0) for m b = at LO but where T(W R ) ^ at NLO. Generally 
speaking, Standard Model radiative corrections can populate density matrix elements 
elements that vanish at leading order. Consequently radiative corrections change angular 
decay distributions as do new physics effects. The lesson to be learned is clear. Before 
ascribing a given polarization effect to new physics one has to make sure that it does 
not result from radiative corrections of old physics. In the above example a non-zero 
F(W R ) 7^ could result from an admixture of a right-handed charged current (new 
physics) or from radiative corrections. Only the precise knowledge of the magnitude of 
the radiative corrections allows one to exclude or conclude for new physics effects. 

In this review we will mostly be concerned with radiative corrections to the current- 
induced transitions t — > b, b — > c and I — > I'. It should be clear that the additional gluons 
(or photons) in radiative correction calculations only couple to the q 1 — » q 2 (or I — > I') 
side of the relevant semileptonic (or leptonic) transitions. This implies that the structure 
of the NLO radiative corrections is the same in the two classes of processes. The current- 
induced transition I — > I' is not in the Standard Model form but in the charge retention 
form. However, the charge retention form can be linked to the Standard Model form by 
the remarkable property that the (V — AY(V — A)^ interaction is an eigenvector under 
Fierz crossing (see Sect. 3). One therefore has 

75Vp7a(]l - 75)^] + h.c. (2) 

lb)v e ][ei a {t - 75)//] + h.c. (3) 

All of the above three transitions are therefore governed by the same matrix elements. 
In particular the NLO QED and QCD radiative corrections have the same structure 
for all three transitions. In Table 1 we provide a list of the processes for which the 
radiative corrections to the rates and the polarization observables can all be calculated 
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Table 1: List of processes based on the same current-induced matrix elements 



parton level particle level references 



ttf)^b + W+a) t(])^X b + W+(]){^l + + v l ) [4,5,6,7] 

t(\)^b + H + t (X)^X b + H+ [8] 

6(T) — c(T) + 1' + vx fl— >X C + £>«(T) [9] 

a 6 (T)^x c + j dW(t) [io] 

6(T) — «(T) + T +n B-^X u + n(p(f)) [9] 

A b (T)^x u + 7r(p(T)) [io] 

l-{])^l'-(]) + vi + vv [3] 
(/,/') = (/i,e), (r,/x), (r,e) 



from essentially the same set of NLO matrix elements. In this review we shall only present 
a few sample results from the processes in Table 1 because the results are to numerous 
to fit into a review of the present size. Instead we attempt to share our insights into the 
general features of spin physics and illustrate these with sample results taken from the 
processes in Table 1. 

2 Miscellaneous Remarks on Polarization Effects 
2.1 Examples of 100% polarization 

In this subsection we discuss examples of 100% polarization. Cases of 100% polarization 
are usually associated with limits when masses or energies of particles become small or 
large compared to a given scale in a process. For example, an on-shell or off-shell gauge 
bosons radiated off massless fermions (leptons or quarks) is purely transverse since fermion 
helicity is conserved by the vector and axial vector Standard Model couplings of the gauge 
bosons. Contrary to this, massive gauge bosons become purely longitudinal in the high 
energy limit (see Sect. 11). We discuss two examples of 100% polarization in more detail. 
Namely, the case of a left-chiral fermion which becomes purely left-handed in the chiral 
limit when to/ — > and the case of soft gluon radiation from heavy quarks. 
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Relativist ic left— chiral fermions 

Consider a fermion moving along the z-axis. The positive and negative helicity polar- 
ization four-vectors of the fermion are given by s± = ±(p; 0, 0, E)/m. At first sight it 
appears to be problematic to take the limit m — > in expressions involving the polariza- 
tion four-vectors because of the denominator mass factor. However, the saving feature is 
that s± becomes increasingly parallel to the momentum when m — > 0. In fact, one 
has 

4 = ±p" + 0(m/£), (4) 

since E = p + (m 2 )/(E + p). Therefore the projectors onto the two helicity states simplify 
to 

1 m->0 1 

u(±)u(±) = -y + m )(l + 75 ^ ± ) ► -^(1 T 75 ). (5) 

This shows that a right /left-chiral fermion is purely right /left-handed in the chiral 
limit. Thus the final state electron emerging from leptonic /x-decay is purely left-handed 
in the chiral limit except for the anomalous spin-flip contribution appearing first at NLO 
which populate also the right-handed state. The anomalous spin-flip contribution will 
be discussed in Sect. 12. 

The approach to the chiral limit is well described by an approximate formula frequently 
discussed in text books. The argument goes as follows. 

Introduce left-chiral fermion spinors according to ul(±) = 1/2(11 — 7s)-u(±), where 
u(±) are helicity A = +1/2, —1/2 spinors. The longitudinal polarization of a left-chiral 
electron is then calculated by taking the ratio of the difference and the sum of the A = 
+ 1/2 and A = —1/2 scalar densities according to 



p i 4(+R(+)-4(-RR f6) 
4(+K(+) + 4(-K(-) 

(E + m-p) 2 -(E + m + p) 2 = _p_ = 

(E + m-p) 2 + (E + m + p) 2 E P ' U 

A corresponding result holds for left-chiral positrons where P\ = f3. 

The accuracy of the above result can be checked with the corresponding expression 
for the Born term polarization of the electron in leptonic /i-decay. One finds (see e.g. 

x(3 -2x + y 2 ) 
x(3 - 2x) - (4 - 3x)y 2 



p i _ x { 6 - ax + y~) 
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where x = 2E e /m tl and y = m e /m^. For the Born term contribution the correction 
to the approximate result (6) P\ = —f3 is of 0(1%) or less such that the correct and 
approximate curves are not discernible at the scale of Fig. 1. Contrary to this the NLO 
polarization deviates substantially from a simple P\ = —(5 behaviour (see Fig. 1) even 
though one has again P l e oc —(3 at NLO [3]. This may have to do with the fact that the 
radiative corrections involve an extra photon emission from the electron leg. The electron 
is therefore no longer left-chiral at NLO as in the LO Born term case. 



Soft gluon or soft photon radiation 

When a soft gluon (or a soft photon) is radiated off a fermion line it is 100% polarized in 
the plane spanned by the fermion and the gluon (or photon). To see this in an exemplary 
way consider the tree graph matrix element of a soft gluon radiated off a top or antitop 
in the process e + e~ — > tig. The soft gluon matrix element reads 

Ji{.s.g.)(i\ _ 71 (Born) [ Pt _ Pt ) * /q\ 

" 1 J M W* Pt-k) a[ } [ ) 

with e*(±) = (0;=FM,0)A/2. 

The 2x2 density matrix of the gluon can be expanded in terms of the unit matrix 1 
and the Pauli matrices <7j according to 



h+h* + h + h*_ \ _ 1 
h-h\ h_h*_ ) ~ 2 



[a-l + Z-a), (10) 



where £ denotes the Stokes "vector" . We have set "vector" in quotation marks because 
the Stokes "vector" does not transform as a vector but rather as a spin 2 object under 
three-dimensional rotations. From (9) one sees that the helicity amplitudes are relatively 
real and that h + = —h-. The normalized Stokes "vector" is then given by 

f/<7= (-1,0,0), (11) 

i.e., in the terminology of classifying the polarization states of a gluon (or a photon), the 
polarization of the gluon is 100% linearly polarized in the production plane. 

In Fig. 2 we show a plot of the linear polarization of the gluon as a function of the 
scaled energy of the gluon for e + e~ — > tig and for e + e~ — > ccg [11]. At the soft end of 
the spectrum the gluon is 100% polarized and then slowly drops to zero polarization at 
the hard end of the spectrum. It is quite remarkable that a high degree of polarization is 
maintained over a large part of the spectrum. 
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Figure 2: Linear polarization of the gluon as a function of the scaled energy of the gluon 
in e+e" -> Q Q g [11]. 

2.2 Examples of zero polarization 

The mean of a single spin polarization observable such as (ff-p) is a parity-odd measure. 
In strong interaction processes, which are parity conserving, there are no single spin 
polarization effects, i.e. (a • p) = (except for the T-odd single spin effects oc (a ■ (pi x 
Pj)) to be discussed later on). This is different in weak and electroweak interactions 
where the presence of parity violations induces a wide variety of single spin polarization 
phenomena. As an example we take e + e~-annihilation into heavy quarks where single spin 
polarization phenomena occur due to (7, Z) -exchange. NLO corrections to these single 
spin polarization phenomena have been discussed in [12, 13, 14, 15, 16, 17]. Another 
example is the semileptonic decay of a polarized heavy quark where the NLO corrections 
to the lepton spectrum have been presented in [18, 19]. 

However, there are double spin polarization effects in strong interactions since (<?•<?) 
is a parity even measure. An example is the correlation between the spins of the top 
and the antitop produced in hadronic collisions [20, 21, 22, 23]. Naturally, double spin 
polarization effects also occur in weak interaction processes such as in high energy e + e~- 
annihilation (see [16, 24, 25, 26]). However, even in weak processes single spin polarization 
effects can vanish when one approaches phase space boundaries. Consider, for example, 
the process e + e~ — > tt at threshold. In general, the vector current induced amplitude is 
s- and d-wave, while the axial vector current amplitude is p-wave. As one approaches 
threshold the s-wave amplitude will dominate and thus there will be no vector-axial 
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vector interference. Therefore single top polarization goes to zero in e + e~ — > ti as one 
approaches the threshold region. 

Another example of zero polarization are NLO T-odd measures in the process e + e~ — > 
qq(g) for massless quarks. T-odd effects arise from the imaginary parts of one-loop 
contributions. However, in the case e + e~ — > qq(g) the imaginary parts can be shown 
to be proportional to the Born term contribution in the mass zero limit. Hence, the 
T-odd measures are zero in this case [27]. When m q ^ the imaginary parts are no 
longer proportional to the Born term and, consequently, one obtains nonzero values for 
the T-odd measures [28, 29]. In the crossed channels (deep inelastic scattering, Drell- 
Yan process) the imaginary parts no longer have Born term structure and one obtains 
non-vanishing T-odd effects even for zero mass quarks [30, 31, 32, 33]. 

2.3 Mass effects 

Mass effects make the radiative correction more complicated. We attempt to illustrate 
the complication brought about by an additional mass scale by listing the (LO + NLO) 
results for the total rate of leptonic \x decays \i~ — > e~ + v e + for (i) y = m e /m ll = 
and for (ii) y = m e /m^ ^ 0. In case (i) one has 




(12) 



For case (ii) one obtains 



r = 



^f((l-^)(l-^ + ^)-12l/ 4 ln(y) + ^{ 
+ 1(1 - y 4 )(75 - 956y 2 + 75y 4 ) - 2y 4 (36 + y 4 ) ln%) 



24 

(l-32y 3 + lQy 4 -32y 5 + y 8 ) 



7T 2 



2 

- ^(60 + 270y 2 - % 4 + 17y 6 )y 2 ln(y) 

- ^(l-y 4 )(17-6% 2 + 17y 4 )ln(l-y 2 ) 

+ 4(1 - y)\l + 4y + Wy 2 + 4y 3 + y 4 ) ln(l - y) ln(y) 
+ 4(1 + y) 4 (l -4y + 10y 2 - Ay 3 + y 4 ) ln(l + y) ln(y) 
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+ 2(3 + 32 2 / 3 + 4% 4 + 32 2 / 5 + 3!/ 8 )Li 2 (-y) 

+ 2(3-32y 3 + 4% 4 -32y 5 + 3y 8 )Li 2 (y)jy (13) 

The length of the NLO term in (13) illustrates but does not describe the added com- 
plication when introducing a new mass scale. The inclusion of full mass effects is not 
of much relevance for leptonic fj,~ — > e~ decays (except in the threshold region) where 
(m e /m^) 2 = 2.34 • 10~ 5 but may be relevant for the leptonic r — > fi decays, where e.g. 
(m^/m T ) 2 = 3.54 • 10~ 3 . In fact, experimental data on leptonic r decays do show a 
mass dependence of the partial rates on the daughter lepton's mass. For the decays 
r~ — > /i~ + + z/ T and r~ — > e~ + z/ e + z/ T one finds branching ratios of 17.37 ±0.06% and 
17.84±0.06% [34]. In agreement with experiment the LO result predicts a reduction of the 
partial rate by 2.82% for the r~ — > yT mode relative to the t~ — > e~ mode. As concerns 
NLO effects, the experimental errors on the two partial rates are still too large to allow 
for checks on the mass dependence of the NLO result. Quite naturally, for semileptonic 
b — > c decays, where y = m c /m b pa 0.3, the mass dependence of the daughter quark must 
be kept. 

Mass effects also populate density matrix elements which are zero when a given mass is 
set to zero. A prominent example is again e + e~ — > qq which is purely transverse for mass 
zero quarks but acquires a longitudinal component when the quark becomes massive. In 
the same vein one has the Callan-Gross relation F 2 = xFi in deep inelastic scattering for 
mass zero quarks which is spoiled when the quarks become massive. Another example 
is the density matrix of the electron in leptonic muon decay. Its off-diagonal elements 
contribute to the transverse component of the electron's polarization (see (20)). The off- 
diagonal elements are proportional to the positive helicity amplitude h±/2 which in turn 
is proportional to the mass of the electron (/11/2 ~ m e ). One therefore has ~ m e [3]. 
A third example is the transition of the top quark into a right-handed W + (Wr) and a 
bottom quark. For m h = the left-chiral bottom quark becomes purely left-handed and, 
from angular momentum conservation, one therefore finds T(Wfi) = (see Fig. 3) (see 
e.g. [5]). 

2.4 Some polarization measures are NLO effects 

Whenever a given polarization component is zero at the LO Born term level, and the 
vanishing of this polarization component is not due to general symmetry principles, it is 
very likely that this polarization component will be populated by gluon or photon emission 
at NLO. We shall list a few examples where this is the case. 
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b 

IXWr) = o 



r(w R ) ± o 



Figure 3: The decay t -> b + W£ at LO and NLO 



A prominent example is the Callan-Gross relation which is violated by gluon emission 
at NLO. Similarly, a longitudinal contribution is generated in e + e~-annihilation into 
massless quarks by gluon emission. T-odd measures derive from imaginary parts of one- 
loop contributions, i.e. they start only at NLO. As mentioned before, a transverse-plus 
polarization of the W + in the decay t — > b + W + is generated by gluon emission as 
illustrated in Fig. 3. At LO the longitudinal polarization of massive quarks in e + e~ — > 
QQ is purely transverse, i.e. the polar distribution of this polarization component is 
proportional to (1 + cos 2 #). The appearance of a longitudinal dependence proportional 
to sin 2 6 is a NLO effect [13]. Anomalous helicity flip contributions involving massless 
quarks or leptons occur only at NLO as discussed in Sect. 12. 

We have already mentioned the T-odd measures which, in a CP invariant theory, 
obtain contributions only from the imaginary parts of NLO one-loop contributions. Apart 
from the examples listed in Sect. 2.2 we mention the case of the transverse polarization 
of a top quark produced in hadronic collisions [35, 36, 37]. 

2.5 Presentation of NLO results of polarization observables 

We have always chosen to present our analytical and numerical NLO results on polariza- 
tion observables in the form 



(P) 



N 



D 



N LO (l + a s N NLO ) 
T> LO (l + a s D NLO ) ' 



(14) 



and not as 



{P) 



N LO 




)) 



(15) 
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An argument for our preference of the form (14) over (15) is the following. It is 
common usage to try and extend a perturbation series beyond a given known order by 
Pade improving the perturbation series. For example, for a second order perturbation 
series the Pade P(l, 1) improvement reads 

\ _|_ Q; "I"" 2 

1 + a s0l + a 2 s a 2 « — (16) 

1 — a s — 

" ai 

What has been done in the Pade P(l, 1) improvement is to fix the first order coefficients 
in the numerator and denominator of the r.h.s. of (16) by expanding the r.h.s. in powers 
of a s and then equating coefficients in (16). Note that at the NLO level Pade's method 
does not yet allow one to reconstruct (14) from (15). Considering the successes of the 
Pade improvement program in other applications we nevertheless believe that the form 
(14) is a better approximation to the whole perturbation series than (15). 

As an example where the use of the form (15) can lead to a misunderstanding is the 
NLO result of the mean forward-backward asymmetry (Afb) i n e + e~-annihilation into a 
pair of massless quarks which, when using the form (15), is sometimes stated as 

LO 

(Afb) = ^£§(1 -a./*). (17) 

This result is very suggestive of a de facto non- vanishing radiative correction to <tfb- 
However, the true result is <Jp^° = [38]. The non-vanishing radiative correction to (Afb) 
in (17) is just a reflection of the well-known NLO result for the rate a = cr LO (l + a s /iv) 
in the denominator of (14). To our knowledge the vanishing of <Jp^° evades a simple 
explanation and must be termed to be a dynamical accident. 



3 Fierz Transformation 

Fierz crossing exchanges two Dirac indices in the contracted product of two strings of 
T-matrices according to (i not summed) 

FPUFmU - [r, M ] a5 [r, M ] 7/3 , (is) 

where the five currents in the set are conventionally labelled by % — S, V, T, A, P with 
r}^ = {1, 7 M , a^, 7^75, 75}. The set of crossed configurations can then be expressed 
in terms of the set of uncrossed configurations. The five-by-five matrix relating the two 
sets is called the Fierz crossing matrix C Fierz- It is clear that one gets back to the original 
configuration when crossing twice, i.e. C Fierz = 1. The eigenvalues A of the Fierz crossing 
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matrix are thus A = ±1 as can easily be seen by going to the diagonal representation. In 
explicit form the Fierz crossing matrix is given by 



C 



Fierz 



1 

4 



V 



1 

4 
12 
-4 

1 



1 

-2 

-2 
-1 



1/2 


-2 


1/2 



-1 
-2 


-2 

1 



1 

-4 
12 
4 
1 



(19) 



/ 



where the rows and columns are labeled in the order % — S, V, T, A, P. 

As (19) shows the trace of the Fierz crossing matrix is Tr(CF ierz ) = —1. This implies 
that the Fierz crossing matrix has three eigenvalues A = — 1 and two eigenvalues A = +1. 
When discussing physics applications one also needs the corresponding Fierz crossing 
matrix for parity odd products of currents. This can easily be derived from multiplying 
the parity even case by 75. 

Amazingly, the charged current interaction of the Standard Model (V — A)^(V — Ay 
is an eigenvector of Fierz crossing, with eigenvalue A = — 1. When relating the Standard 
Model and charge retention forms of the Lagrangian as in (2) the minus sign from Fierz 
crossing is cancelled from having to commute the Fermion fields an odd number of times 
in order to relate the two forms. 

In pre-QCD days it was considered unfortunate that the eigenvalue of the (V— A)^— 
A)^ current-current product is negative. If it were positive one would have a very natural 
explanation of the so called A J = 1/2- or octet rule in weak nonleptonic decays. To see 
this consider the direct products 3® 3 = 3 a © 6 S and 3 © 3 = 3 a © 6 S . With the wrong 
eigenvalue A = +1 (and the minus sign from Fermi statistics) one would then remain with 
3 a © 3 a = 1 © 8 and one thus would have explained the famous octet rule for nonleptonic 
weak interactions. Nevertheless, with the advent of QCD and colour, it was realized that 
an important class of diagrams in nonleptonic baryon transitions between ground state 
baryons was subject to the octet rule due to the simple Fierz property of the nonleptonic 
current-current product. This discovery is sometimes referred to as the Pati-Woo theorem 
[39] even though the discovery of Pati and Woo was predated by the paper [40]. 



4 Counting Spin Observables 

When setting up a problem involving the spin of particles it is always instructive to 
first denumerate the complexity of the problem and count the number of independent 
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Table 2: Independent components of the single spin density matrix 



spin 



J 



rate real imaginary sum 
trace J(2 J + 3) J(2 J + 1) (2 J + l) 2 



J = 1/2 

J= 1 
J = 3/2 



1 



1 



1 



2 



9 



5 



3 



6 



1 



16 



4 
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structures of the problem. A particularly efficient and physical way to do so is to count 
the number of independent elements of the spin density matrices involved in the process. 

The single spin density matrix p\y is a (2J+1)(2J+1) hermitian matrix (pw = Py\)- 
It thus has (2 J + l) 2 independent components of which (J + 1)(2J + 1) are real and 
J (2 J + 1) are imaginary. The rate is represented by the trace of the density matrix. In 
Table 2 we list the corresponding degrees of freedom for the first few spin cases J = 1/2, 1, 
3/2. The density matrix is also frequently represented in terms of its angular momentum 
content. For the three cases listed in Table 2 the angular momentum content is given by 
1/2(8)1/2 = 0©1, 1©1 = 0©1©2 and 3/2©3/2 = 0©1©2©3 where the one-dimensional 
representation "0" is equivalent to the trace or rate. It is clear that one obtains the correct 
respective number of spin degrees of freedom (2 J + l) 2 in Table 2 when adding up the 
dimensions of the angular momentum spaces appearing in the decomposition J © J = 
0©... ©2J. 

In the spin 1/2 case the four spin degrees of freedom are the rate, the two real com- 
ponents and the imaginary component of the spin 1/2 density matrix. The unnormalized 
density matrix p is usually written as 




(20) 



In the helicity basis the four spin degrees of freedom are labelled by 



a 



rate , 
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ii. 



if 



ii. 



•perpendicular" polarization P , 
normal" polarization P N , 



longitudinal" polarization P . 



(21) 



The y-component P can be seen to be a T-odd measure since one needs a vector product 
to construct the normal to a given plane. For example, in /z-decay P N is determined by 
the expectation value (a e ■ (<t m x p e )) which is a T-odd measure. In the context of this 
review P N obtains contributions only from the imaginary parts of one-loop diagrams. For 
all processes listed in Table 1 in Sect. 1 the one-loop diagrams are real and thus P N = 
in these processes. This is different in e + e~ — > qqg where the one-loop diagrams possess 
imaginary parts and where one therefore has P N ^ (see e.g. [14, 41]). 

Returning to Table 2 we shall see that not all (2J + l) 2 spin degrees of freedom are 
accessible (angular momentum conservation) or measurable (parity) in general. 

In the following we will be concerned with double spin density matrices. Their spin 
degrees of freedom are given by the products of the single spin degrees of freedom. For the 
cases discussed in this review the double density matrix is a sparsely populated matrix, 
due to angular momentum conservation and the absence of imaginary parts from loop 
contributions. Take, for example, the double spin density matrix of the decay t(|) — > Xf,+ 
W + (f). On naive counting one would expect N = (1 (trace) + 20(real) + 15(imaginary)) 
spin degrees of freedom. This number is considerably reduced by considering the following 
two facts. 

1. The one-loop amplitude does not possess an imaginary (absorptive) part. This can 
be seen by taking a look at the one-loop vertex corrections for the transitions dis- 
cussed in Table 1. The vertex correction does not admit of real on-shell intermediate 
states, i.e. it does not have an absorptive part. The 15 imaginary components of 
the double density matrix are zero. 

2. Angular momentum conservation. In the rest frame of the top the decay into X5 
and W + is back-to-back and thus anti-collinear. Since one is summing over the 
helicities of the X b one has \x b = \' Xb and thus A t — A^ — X' t — X' w (see Fig. 3). 
Taking this constraint into consideration one remains with six diagonal and two 
non-diagonal double spin density matrix elements (see Table 3) where one refers 
to diagonal and non-diagonal elements when (A t , \w) = (AJ, A'^) and (A t , \w) 7^ 
(\' t , \' w ), respectively. One has non-diagonal transitions for the two cases A t — \w = 
±1/2. Concerning the non-diagonal density matrix elements it is easy to see from 
the master formula (26) to be discussed in Sect. 5 that they generate azimuthal 
dependences in angular decay distributions. 
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We shall now go through the exercise and count the spin degrees of freedom by a differ- 
ent method, namely by counting the number of covariants of the process. We emphasize 
that counting spin degrees of freedom via density matrix elements is generally safer than 
counting the number of independent covariants because there exist nontrivial identities 
between covariants in D = 4 dimensions. The decay — > X b + W + ('\) will illustrate 
this point. 

Consider the expansion of the spin dependent hadron tensor of the process into co- 
variants constructed from the metric tensor, the Levi-Civita tensor, the independent 
momenta of the process and the spin four-vector of the top. One has 

H v> = (-g'v'^+rfplHz-ie^pt^Ha) 

- (q-s t )(- <T d + ftft G2 - ie^pt^a G 3 ) 

+ (sW t + sffi) G 6 + ie^p tp s ta G 8 + ie^q p s ta G 9 . (22) 

There are nine covariants and associated with it nine invariants which obviously does 
not agree with the number eight counted before using helicity counting. The discrepancy 
arises because there exist a nontrivial four-dimensional identity between three of the nine 
covariants which reads 

q-s t ^ vpa Pt, P qa - q 2 ^ pa Pt, P s t a + q-Pt e^q p s t ^ = 0. (23) 
This identity can be derived using the Schouten identity which reads 

^/i[/iiM2/*3M4M5] = 0' (24) 
where the symbol "[. . .]" denotes antisymmetrization and where 

+ cycl.(/ii, ^2,^3, ^4,^5). (25) 

The Schouten identity is just the statement that it is impossible to place four index 
values in an antisymmetric fifth rank tensor, or in the language of Young Tableaux, that 
a Young Tableau with five vertical boxes is identically zero in four dimensions. 

This illustration is not only of academic interest but has numerical implications since 
one can get into a terrible mess numerically if one works with a redundant set of covariants 
and tries to do matrix inversions involving the overcounted set of degrees of freedom. 
Although the Schouten identity seems rather obvious nowadays there have been examples 
in the literature where Schouten-type of identities have been overlooked. 
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Table 3: Non-vanishing double density matrix elements in t(f ) — > Xf, + W^ + (f) 



A* 


Xw 


Xt — Xw 


1/2 


1 


-1/2 


1/2 





1/2 


1/2 


-1 


3/2 


-1/2 


1 


-3/2 


-1/2 





-1/2 


-1/2 


-1 


1/2 



We conclude this section by enumerating the number of spin degrees of freedom in 
the various decay processes discussed in this review. The processes are listed in Table 4 
together with the number of spin degrees of freedom. The corresponding references to the 
papers in which these decays were treated can be found in Table 1. In this list we have 
taken into account that the decay D* — > 71(7) is parity conserving. In parantheses we list 
the number of T-odd observables for each process which are zero due to the absence of 
imaginary contributions in the one-loop vertex correction. 

5 Angular Decay Distributions 

It should be clear that one needs to do polarization measurements in order to disentangle 
the full structure of particle interactions. Polarization measurements are particularly 
simple when the particle whose polarization one wants to measure decays. The angular 
decay distribution of the decay products reveals information on the state of polarization 
of the decaying particle. The information contained in the angular decay distribution is 
maximal when the particle decay is weak. The fact that the angular decay distribution 
reveals information on the polarization of the decaying particle is sometimes referred to 
as that the particle decay is self-analyzing. 

There are principally two ways to obtain angular decay distributions which we will re- 
fer to as the non-covariant and the covariant methods. In the non-covariant method one 
makes use of rotation matrices whereas in the covariant method one evaluates scalar prod- 
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Table 4: Number of measurable double density matrix elements for the processes listed 
in Table 1 



process 

t(T) — > 6 + T^+(t) 

t(T) ^b + H+ 
B --^X c + D s 
B ^ Xc + D* s (t) 
A 6 (T) — > X c + P s 

a 6 (t)^x c + j d:(t) 
i(t) — m) + ^ + ^ 



spin degrees of freedom 

AT = 1 + 7 (+2) 
AT = 1 + 1 
A" = 1 
N =1 + 1 
N =1 + 1 
AT = 1 + 4 (+1) 
AT = 1 + 4 (+1) 



ucts of four-vectors involving momenta and spin four-vectors in given reference frames. 
We shall discuss one example each of the two methods. 

As an example for the non-covariant method we write down the angular decay distri- 
bution of polarized top decay — > b + W + followed by W + — > l + + v\. The angular 
decay distribution can be obtained from the master formula 



where pa 4 a{(^p) is the density matrix of the top quark which reads 

1 / 1 + P cos Op P sin e P \ 
2\ Psin0 P 1-PcosOp) 

P is the magnitude of the polarization of the top quark. The H**^, are helicity matrix 
elements of the hadronic structure function H^. The sum in (26) extends over all values 
of Xw,X' w ,X t and X' t compatible with the constraint Xw — X' w = X t — X' t . The second 
lower index in the small Wigner d(#)-function d\ wl is fixed at m — 1 for zero mass 
leptons because the total m-quantum number of the lepton pair along the / + direction is 
m — 1. Because there exist different conventions for Wigner's d-functions we explicate 
the requisite components that enter (26): d^ = (1 + cos#)/2, dj x = sm9/\^2 and d 1 _ 11 = 
(1 — cos#)/2. The fact that one has to specify the phase convention of the Wigner's 
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(^-functions points to one of the weaknesses of the non-covariant method to obtain the 
correct angular decay distributions: one has to use one set of consistent phase conventions 
to obtain the correct signs for the angular factors. For someone not not so familiar with the 
angular momentum apparatus this is not always simple. For the polar angle dependencies 
the correctness of a sign can always be checked by using physics arguments. This is more 
difficult for the azimuthal signs. 

Including the appropriate normalization factor the four-fold decay distribution is given 
by [4, 6] 



dT 1 G F \V tb \ 2 m 2 

dq d cos 9 pd cos 9d<j) 4ir \f2ii 
3 



w 



+^(Hu + Pcos9 p H uP )(l + cos 2 9) 

8 

+ ^{H L + P cos 6 p H L p ) sin 2 9 + - (H F + P cos 6 p H F p ) cos 6 
3 „ . „ „ 3 

4 



3 1 
r P sin 9 p Hjp sin 26 cos H — -j=P sin 9 p H A p sin 9 cos > , (28) 

2y2 v2 J 

where the helicity structure functions Hy, H F etc. are linear combinations of the helicity 
matrix elements HJ We have taken the freedom to normalize the differential rate 
such that one obtains the total t — > b + W + rate upon integration and not the total rate 
multiplied by the branching ratio of the respective W + decay channel. 

The polar angles 9p and 9, and the azimuthal angle that arise in the full cascade- 
type description of the two-stage decay process — * W + {-^ l + + +X b are defined in 
Fig. 4. For better visibility we have oriented the lepton plane with a negative azimuthal 
angle relative to the hadron plane. For the hadronic decays of the W into a pair of light 
quarks one has to replace by (q, q) in Fig. 4. We mention that we have checked 

the signs of the angular decay distribution (28) using covariant techniques. 

At first sight it seems rather strange that the angular analysis is done in two different 
coordinate systems, namely the top quark rest system and the W rest system. This runs 
counter to common wisdom that invariants should always be evaluated in one reference 
system. Some insight into the problem may be gained by using the orthonormality and 
completeness relation of polarization vectors to rewrite the contraction of the hadron and 
lepton tensors in a form which exhibits the correctness of using two different reference 
systems for the cascade decay. The orthonormality and completeness relation read 



Orthonormality: V e *''( m ) eV ( m ') = W' m ? m' — S, ±,0, (29) 
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Figure 4: Definition of the angles Op, and in the cascade decay — > + and 
-> /+ + i/, 

Completeness: ^<?(m)e* v {m')g mm > = ^ Ji, v = 0, 1, 2, 3, (30) 

m,m'=S,±,0 

where sw = diag(+, -, -, -), (m,m' = S,±,0) and g» v = diag(+, -, -, -), = 
0, 1, 2, 3). The scalar or time component of the polarization four-vector is denoted by S. 
On using the completeness relation one then has 

ITH^ = L^'^H^ = L^ v ,6»\m)6*»{m')g mm ,6* v \ny{n')g nn ,H^ 

= (L^\m)e^\n))(H^(mY(n'))grn m '9nn'. (31) 

The point is that the two Lorentz contractions appearing on the second line of (31) can be 
evaluated in two different Lorentz frames. The leptonic invariant L^yit^ (m)e* u (n) can 
be evaluated in the (lu) CM frame (or in the W + rest frame) while the hadronic invariant 
H^ l/ e* fl (m f )e u (n f ) can be evaluated in the rest frame of the top quark. Another advantage 
of this method is that one can easily incorporate lepton mass effects as is mandatory for 
b — > c or the rare b — > s transitions if the charged lepton or leptons in the final state are r 
leptons. This technique was used to derive angular decay distributions including leptonic 
polarization effects for semileptonic and rare bottom meson decays in [42] and [43]. Note 
also that the correct phase choice for the polarization vectors is no longer crucial since 
the polarization vectors always appear as squares in the orthonormality and completeness 
relations. 
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As an example of the covariant method we discuss the angular decay distribution of 
polarized fi decay into a polarized e. From helicity counting as described in Sect. 4 one 
knows that there are altogether five spin-dependent structure functions and one spin- 
independent structure function describing the leptonic decay of a polarized muon into 
a polarized electron. We thus define a spin-dependent differential rate in terms of six 
invariant structure functions A i . Accordingly one has [3] 



+ A 5 (s^s e ) + -^A 6 e af3j sPZp P eSls 5 e . (32) 

m ix 

Eq. (32) will be evaluated in the rest system of the muon where p^ = (m^; 0,0,0) 
and p e = (_E e ;0,0, \p e \) = (m )l /2)(x; 0, 0, x(3). The velocity of the electron is denoted by 

f3 = \J\ — 4y 2 /x 2 where y = m e /m^ and x = 2E e /m ll denotes the scaled energy of the 
electron. In the rest frame of the \i~ the polarization four-vectors of the \i~ and are 
given by 



s 



(0;CJ, (33) 



a ,^e'Pe -> ^e'Pe -> \ /oa\ 

S e = ( ;™e+ 7TT— tPe), (34) 

m e m e (E e + m e ) 

where the polarization three-vector of the fx~ and the quantization axis n e of the spin 
of the e~ in their respective rest frames read (see Fig. 5) 

C„ = (sni0p,O,cos0p) (35) 

and 



n e = (sin 9 cos x, sin $ sin x, cos 6) . (36) 

Eq. (35) holds for 100% polarized muons. For partially polarized muons with mag- 
nitude of polarization P the representation (35) has to be multiplied by P such that 

P, = PC,- 

The scalar products in (32) can then be evaluated with the result 
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Figure 5: Definition of the angles 9p, 9 and x m polarized muon decay 



Pe' s u = 7T X (3P COS 6p, 





(37) 



Including the correct normalization one finally arrives at the angular decay distribution 



where the relation between the invariant structure functions Ai and the frame-dependent 
spectrum functions G{ can be found in [3]. G\ is the unpolarized spectrum function, 
G 2 and G 3 are single spin polarized spectrum functions referring to the spins of the 
/i~ and e~, resp., and G 4 , G 5 and G 6 describe spin-spin correlations between the spin 
vectors of the muon and electron. G e represents a so-called T-odd observable. This is 



[3] 



dT 



fix r (Gi + G 2 P cos 9 P + G 3 cos 9 + G 4 P cos 9 P cos 9 



dx dcos 9p 



+ G 5 P sin 9 P sin 9 cos \ + GePsin^psin^sinx), 



(38) 
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evident when rewriting the angular factor multiplying G e in (38) in triple-product form, 
i.e. sin dp sin 9 sin x = |Pe| _1 Pe • (Cm x ^e)- As stated before Gq is identically zero in the 
Standard Model since, on the one hand, the weak coupling constant Gf is real and, on the 
other hand, the loop contributions do not generate imaginary parts. In the terminology of 
Sect. 4 G\— G 4 are diagonal structure functions, while the azimuthally dependent structure 
functions G 5 and Gq are non-diagonal. As mentioned before, the structure function G 5 
is proportional to the electron mass and therefore vanishes for m e — > 0. Naturally, the 
angular decay distribution (38) can also be derived in the helicity formalism. 

In the two examples discussed in this section we chose to orient the z-axis along the 
three-momentum of one of the particles. In this case one says that the angular analysis 
is done in the helicity system. Of course, other choices of the orientation of the z-axis 
are possible. A popular choice is the transversity system where the z-axis is normal to a 
given plane spanned by the momenta of the decay products. A more detailed discussion 
of the choice of frames can be found in [42] where the angular decay distribution in the 
semileptonic decays B — > (D, D*) + I + v including lepton mass effects was investigated. 
We mention that the transversity system is the preferred choice when discussing the CP 
properties of the final state [44]. 

6 One-Loop Amplitude 

We present our results in terms of the three vector current amplitudes Fj [i = 1,2,3) 
and the three axial vector current amplitudes Ff- [i = 1,2,3) defined by ( Jj = qbl^Qt, 

(b(pb)\tf\t(pt)) = u b (Pb){l^ +pt, li F? + Pb ,^}u t (p t ), (39a) 

(b(pb)\tf\t(pt)) = uM^F^+p^+p^F^MPt)- (39b) 

In this Section we choose to label the current transition according to the transition t — > b 
as in [6]. As emphasized before, up to the colour factor identical expressions are obtained 
for the /j, — > e transition if the leptonic four-fermion interaction is written in the charge 
retention form. For the form factors one obtains [6] 

2Li 2 (l -wl) -2Li 2 (l - — ) (40a) 



FX 



mt + mi 
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A 4 



2 2 



>±-< In - 4 + to^,) - ( "" + ' 2 "*>!-^ 

2g 2 Vm t V p V 2g 2 m\\f\ 

v = a. m, - m 6 r _ /m t + 2m 6 _ m? - m 2 X ^ /rn| 
47T g 2 I \m t — m b q 2 J Km* . 

\H\ m b q 2 + (m t - m 6 ) (3m t + m b ) 



m 2 - 



ln(w!W M )|, 



g 2 m t - m 6 m 2 ^ 
F 3 y = F^{m u m b ) = F?(m b ,m t ), (40c) 

where we have denoted the scaled (small) gluon mass by A = m g /m t . We define A = 
1 + x A + y A — 2x 2 y 2 — 2x 2 — 2y 2 (x = q 2 /m 2 , y = m b /m t ) and use the abbreviations 

x l — x 2 + y 2 — y r \ x l — x 2 + y 2 — \f\ 

Wi = =, w u = =. 41 

y l + x 2 -y 2 + V\ M y l + x 2 -y 2 -^f\ 

The axial vector amplitudes can be obtained from the vector amplitudes by the 
replacement m t — > — m t) i.e. one has F t A (m t ) = F^(—m t ) (i = 1,2,3). Our one-loop 
amplitudes are linearly related to the one-loop amplitudes given in [45] after correcting 
for a typo in [45] (see also [46]). 

Note that the infrared singularities proportional to In A and the would-be mass singu- 
larities (also called collinear singularities) proportional to In m b all reside in the Born term 
form factors F^ and Ff. They are eventually cancelled by the corresponding singularities 
in the tree graph contribution. A look at the arguments of the log and dilog functions in 
40a shows that the one-loop contribution is purely real as remarked on earlier. 

In the case of the electroweak radiative corrections to t — > b + W + there are altogether 
18 different vertex correction diagrams in the Feynman-'t Hooft gauge [7] as compared 
to the one vertex correction diagram in the QCD and the //-decay cases discussed in this 
section. In addition to the massive one-loop three-point functions one has to calculate 
the many massive one-loop two-point functions needed in the renormalization program. 
We have recalculated all one-loop contributions analytically and have checked them ana- 
lytically and numerically with the help of a XLOOPS/GiNaC package that automatically 
calculates one-loop two-point and three-point functions [47] . Our one-loop results agree 
with the results of [48]. The results are too lengthy to be reproduced here in analytical 
form. The full analytical results will be given in a forthcoming publication [49]. 
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7 Tree-Graph Contribution 



As emphasized earlier on, the NLO tree graph contributions for the QED and QCD cases 
are identical up to NLO except for a trivial colour factor in the case of QCD. Differences 
set in only at NNLO where three-gluon coupling contributions come in. We choose to 
present and discuss the NLO tree graph contributions in the QED radiative corrections 
to leptonic /x-decay as written down in [3]. 

We begin with the Born term contribution. This is an exercise that everyone has 
probably gone through before in the unpolarized case. The point that not everyone is 
familiar with is that it is very simple to include the spin of fermions at the Born term 
level. All that is needed is the substitution p — > p for the fermion's momenta where the 
notation p is explained in the following. The p —> p rule is best explained by looking at 
the trace expression of the charge-side tensor which reads 

Ctm = ^Tr{(y e + m e )(l + l5 ^h a (l - 75 )(^ + m,){t + ^^{t - 75)}. (42) 

The dependence on the polarization four-vectors of the n~ and e~ has been retained in 
(42). 

Since only even-numbered 7-matrix strings survive between the two (11 — 7 5 )-factors 
in (42) one can compactly write the result of the trace evaluation as 

Ctm = Kf.Pe +P»P P e- <f* ' PfPe + 16^%^), (43) 

where 



Pi = P a ,-m»s«, (44a) 
P a e = P a e -m e s a e , (44b) 

and where and s" are the polarization four-vectors of the fi~ and e~. 

In the QED case one has the simplifying feature that the contribution of the last 
term antisymmetric in (a(3) can be dropped. The reason is that the dependence on the 
momentum directions of the V e - and z/^-neutrinos is completely integrated out in the 
differential rate. Thus the neutrino-side of the interaction can only depend on the spatial 
piece of the second rank tensor build from the momentum transfer to the neutrinos (the 
neutrinos are treated as massless) which is symmetric in (aft). Upon contraction with the 
charge-side tensor the antisymmetric piece drops out. 
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Unfortunately the p — > p trick no longer works at NLO order. But still, by using the 
p-notation, the result can be presented in a very compact form. The result for the NLO 
p, — > e tree graph contribution reads 



C<# = V M a M^- e2 \( k ' pe ~ m2e I P»-PA ka Pv + k/3 P»- k -P^ 



7— spin 



+ 



k-Pe k-Pf, J k-p e 

k -Pm + ml p e -p^\ k a f e + fc^ - k-p e g a/3 



k-Pf, k-p e J k-Pf, 



(k-p e )(k-p lx ) ■ ' (k-p e )(k-p^ 

( , - N + Pf P^ - Pe-P^g a(3 , _ . P^fe + P^P" ~ P^PeQ^ 

+ <**> (FrF ( *' p " ) i^F 

The momentum of the radiated photon is denoted by fc. For the aforementioned reason we 
have dropped the contribution of antisymmetric terms. When the antisymmetric e-tensor 
pieces are kept, and when one replaces e — > g s and 1 — > N c Cp = 4 to account for colour, 
one recovers the NLO QCD corrected hadronic tensor for the t —> b transition listed in 

[6]- ' 

In the last line of (45) we have isolated the infrared singular piece of the charge- 
side tensor which is given by the usual soft photon factor multiplying the Born term 
contribution. Technically this is done by writing 

C (a)a/3 = (c {a)af3 - C (Q)Q/3 (softphoton)) + C (a)a/3 (softphoton). (46) 

The remaining part of the charge-side tensor in (45) is referred to as the hard photon 
contribution. It is infrared finite and can thus be integrated without a regulator photon 
mass. 

Analytic phase space integrations without a regulator photon mass are much simpler. 
The regulator photon mass would introduce a new mass scale into the problem which, as 
emphasized before, would complicate the phase space integrations. 

In the phase space integration over the photon momentum the infrared singular piece 
is regularized by introducing a (small) photon mass resulting in a logarithmic mass diver- 
gence in the photon mass. Since the integrand of the soft photon piece is much simpler 
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the phase space integration can be done analytically by carefully taking the appropriate 
mass zero limits in the integrations. The infrared divergence shows up as a logarithmic 
mass divergence in the photon mass. In addition, since the soft photon piece factors the 
Born term tensor, the Born term tensor can be pulled out of the phase space integrations. 
The resulting singular soft photon piece is therefore universal in the sense that it is the 
same for all spin structure functions, i.e. once the soft photon integration has been done 
for the rate the work is done. Eventually the infrared singular piece is cancelled by the 
corresponding singular piece in the one-loop contributions. 

The approximation where only the soft photon (or gluon) piece is retained in the 
NLO tree-graph hadron tensor is called the soft photon (or gluon) approximation. From 
what has been said before it is clear that a NLO radiative correction calculation is much 
simplified in the soft photon (or gluon) approximation. In the next section we shall take 
a specific example, namely the process e + e~ — > tt(g), to investigate the quality of the 
soft gluon approximation. 




Figure 6: Tree-level Feynman graphs for t — > b + W + + 7 in the Feynman-'t Hooft gauge 

It is interesting to note that the NLO hadronic tensor in the elektroweak corrections 
to the process t — > b + W + is only marginally more complicated than the corresponding 
tree graph contribution (45) although the number of contributing diagrams has doubled 
to four [7] (see Fig. 6). We do not list the finite piece of the hadronic tensor in this 
review but only write down the soft photon factor which multiplies the Born term tensor 
just as in (45). When discussing the decay t — > b + W + one has to of course include the 
antisymmetric e-tensor piece in the Born term tensor. The soft photon contribution reads 



_gMi^(s-ph.)^(s.ph.) 
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(47) 
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where the soft photon amplitude is given by (see e.g. [51]) 



" \ P fk Pb -k q -k r»- 1 8J 

Qt = 2/3, <5b = —1/3 and Q w = 1 are the electric charges of the top quark, the bottom 
quark and the VT-boson, resp., in units of the elementary charge e. Gauge invariance of 
the soft-photon amplitude is easily verified when replacing e* — > k^ in (48). It is then just 
a statement about charge conservation Qt~Qb~Qw — 0. It is noteworthy that by setting 
Qt = Qb = 1 an d Q w = in (47), one recovers the soft photon contribution in (45). In 
fact, the whole NLO electroweak tree-graph tensor in [7] reduces to the corresponding 
NLO QED tensor with the above charge replacements. In the same vein, the replacements 
e — > g s and 1 — > N c Cf = 4, to account for colour, and the replacements of the charge 
factors by Q t = Qb = 1 and Q w = will bring one from the electroweak case to the QCD 
case. 

8 NLO Radiative Corrections to e + e~ — ► ti(g) in 
the Soft Gluon Approximation 

In the soft gluon (or photon) approximation one keeps only the soft gluon (or photon) 
piece in the tree graph contribution but includes the full one-loop contribution. New 
structure is thus only generated by the one-loop contribution since the soft gluon piece 
has Born term structure. The structure contained in the hard part of the tree graph 
contribution is lost in the soft gluon approximation. A brief glance at the corresponding 
hard photon piece in (45) shows that it is considerably more difficult to do the analytic 
phase space integration for the hard part than for the soft photon piece. In particular, 
the integration of the hard photon piece has to be done separately for each density matrix 
element. Contrary to this, the integration of the soft photon piece has to be done only 
once since the Born term factor can be factored out of the integral. Technically, the 
soft gluon (or photon) approximation is much simpler than the full calculation if done 
analytically. Of course, if the calculation is done numerically, the integration of the hard 
gluon (or photon) part causes no additional problems. When the soft gluon (or photon) 
approximation is used this is done at the cost of loosing interesting structure contained 
in the hard gluon (or photon) part. 

In order to be able to judge the quality of the soft gluon approximation we have 
calculated the rate for e + e~ — > tt(g) and have compared the results of the full calculation 
with the results of the soft gluon approximation using different cut-off values for the gluon 
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Figure 7: Total cross section for e + e~ — > tt(g) as a function of the scaled gluon energy 
cut \/\max for different cm. energies. Full line: full calculation; dashed line: soft gluon 
approximation 

energy [52]. Fig. 7 shows that the soft gluon result is below the full result for the whole 
range of cut-off values A = E g /y/q^ (\ m ax — 1 — 2m t / A/tp). For maximal cut-off values 
A = \max the soft gluon approximation is 3.4%, 5% and 11.5% below the full result at 
400 GeV, 500 GeV and 1000 GeV, respectively. As emphasized above, much of the rate 
that is being missed by the soft gluon approximation has interesting structure. 

9 Unpolarized Top Decay t — ► b + W + 

We have already discussed various aspects of the decay t — > b + W + in previous sections. 
In this section we concentrate on unpolarized top decay. In particular, we want to discuss 
the mass dependence of the longitudinal piece of the W + boson the measurement of which 
could lead to an independent determination of the mass of the top quark. 

The angular decay distribution for unpolarized top decay can be obtained by setting 
P = in (28). The contribution of the three remaining structure functions Hi, Hjj and 
Hp can be disentangled by a measurement of the shape of the lepton energy spectrum from 
top decay. Given enough data, one can hope to determine the longitudinal contribution 
with 1% accuracy [50]. 
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At LO the mass dependence of the longitudinal contribution is given by T^/T = 
\/(\+2{m w /m t ) 2 ) which gives T L /T = 0.703 using m t = 175 GeV and = 80.419 GeV. 
NLO corrections to the longitudinal contribution were calculated in [5, 6] (QCD) and in 
[7] (electroweak and finite width). Curiously enough the electroweak and finite width 
corrections tend to cancel each other in the structure functions. In Fig. 8 we show the top 
mass dependence of the ratio T L /T. The Born term and the corrected curves are prac- 
tically straight line curves. The horizontal displacement of the two curves is ~ 3.5 GeV. 
One would thus make the corresponding mistake in the top mass determination from a 
measurement of Tl/T if the Born term curve were used instead of the corrected curve. If 
we take m t = 175 GeV as central value, a 1% relative error on the measurement of T L /T 
would allow one to determine the top quark mass with « 3 GeV accuracy. Such a top 
mass measurement would be a welcome alternative to the usual invariant mass determi- 
nation of the top quark mass since the Tl/T measurement is a completely independent 
measurement of the top quark mass. 
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Figure 8: Top mass dependence of the ratio Tl/T (full line: LO, dashed line: correc- 
tions including QCD, electroweak (GV-scheme), finite-width and (mj ^ 0) Born term 
corrections) 
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10 The Decay t(T) -> b + H+ 



Certain extensions of the minimal Standard Model with two Higgs-doublets (as e.g. the 
Minimal Supersymmetric Standard Model) contain charged Higgs bosons such that the 
decay i(t) — > b + H + would be possible if kinematically allowed. The Lagrangian for the 
decay reads 

C = V b (x)(a + lry 5 )%(x)(f> H +(x) + h.c. . (49) 
For the differential decay rate one obtains 

= i (r + pr? ms e P ) = ir(i + p a[ , ms e P ), (to) 

where dp is the angle between the polarization vector of the top quark and the Higgs. 
The cos#p-dependence in (50) is determined by the asymmetry parameter an = r/Tp. 
At LO one obtains {m h = 0) 

®h = 2 , 2 - (51) 

At this point we want to emphasize that one needs the polarization information con- 
tained in (50) to be able to determine the relative sign of the two coupling constants a 
and b. The rate is proportional to (a 2 + b 2 ) and is therefore insensitive to the relative sign 
of the coupling constants. 

We now specify to the so called model 2 where the coupling structure is expressed in 
terms of tan/3 = v^/vi, and where V\ and v-i are the vacuum expectation values of the 
two neutral components of the two Higgs doublets. In model 2 the coupling constants are 
given by 

V t b(m t cot /3 + nib tan/3), (52) 



2^/2m w 

b = — j^- — Vtb{m t cot/3 — m;,tan/3). (53) 
2 v 2m w 

The weak coupling factor g w is related to the usual Fermi coupling constant Gp by g w = 
c lm w <f2^fG F ~. In Fig. 9 we show a plot of the asymmetry parameter olh as a function of 
the mass ratio m^/mt for a fixed value of tan/3 = 10. The asymmetry parameter is large 
and negative over most of the range of mass ratios. The radiative corrections can be seen 
to be quite substantial. 
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Figure 9: Asymmetry parameter an for model 2 with mj = 4.8 GeV and m t = 175 GeV 
(LO: full line, NLO: dashed line) as function of m#+/ m i and tan/3 = 10. The barely 
visible dotted lines show the corresponding m& — > curves. 

11 Goldstone Equivalence Theorem 

Consider the decay t — > 6 + J^(q 2 ) where the current J^(q 2 ) can be thought of as repre- 
senting an off-shell containing spin 1 and spin pieces. Taking g M = (q ; 0, 0, |g|) the 
longitudinal polarization vector of the spin 1 piece is given by e M (L) = (\q\; 0, 0, qoj/y/q 1 
while the scalar polarization vector reads e' 1 (S) = q^/y/q 1 . In the limit \fq^jm t the 
longitudinal polarization vector becomes increasingly parallel to its momentum. In fact, 
one finds e M (L) = q^/^/q 2 + O^^/q 2 / 'q ). In this limit the amplitude A(L) dominates over 
the transverse amplitudes A(T) since A(T)/A(L) ~ \fq I /q$- At the same time one finds 
A(S) ~ A(L) in this limit. This provides a useful check on the high energy limits of the 
relevant unpolarized and polarized rate functions. 

Next consider the scalar projection of the t — > b current transition at lowest order. 
Using the Dirac equation one obtains (q — p t — Pb) 

q^u^il - 75)144 = (p t - p h )^u h ^{t - 75)1*4 

= u b (m t (l + ^ 5 ) +m b u t (l- 75))^*. (54) 

The second line of (54) has the coupling structure of a Goldstone boson in the Standard 
Model. It is a very instructive exercise to follow this argument through to one-loop 
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order including renormalization [53]. In fact, this generalizes to any order in perturbation 
theory: the scalar piece of the charged transition current has the coupling structure of a 
Goldstone boson. Together with the above statement that the longitudinal rate dominates 
at high energies and that it is equal to the scalar rate in this very limit one arrives at 
the Goldstone Boson Equivalence Theorem. In the high energy limit the coupling of a 
gauge vector boson to a fermion current is equivalent to that of a Goldstone boson. This 
simplifying feature has frequently been used in the literature to dramatically reduce the 
effort needed in the computation of high energy processes involving vector gauge bosons. 
For example, the Goldstone equivalence theorem was used in [54] to calculate the bosonic 
two-loop electroweak radiative corrections to the decay H — > 7 + 7 in the limit of a 
large Higgs mass. The same approximation was used to calculate the dominant bosonic 
two-loop electroweak radiative corrections to the decay of a heavy Higgs into pairs of W 
and Z gauge bosons [55]. 

As concerns the coupling of charged Higgs bosons in t — > b + H + discussed in Sect. 10 
the coupling structure of the charged Higgs is equal to that of a Goldstone boson in the 
Standard Model if one specifies to the so called model 1 with cot/3 = 1 (see e.g. [8]). In 
this case one has 

r t ^ b+H + ~ T t ^ b+W + as m w +,m H +/m t — > 0. (55) 

The same relation holds true for the corresponding polarized top rates. All of the 
statements made in this section have been explicitly verified using our NLO results on 
t -> b + W+ [6] and t -> b + H + [8]. 

12 Leptonic Decays of the /x and the r and Anoma- 
lous Helicity Flip Contributions 

We have already discussed various aspects of the leptonic decays of the \x and the r lepton 
in previous sections. A complete solution to the problem has been given in [3] where, for 
the first time, all mass and polarization effects have been included at NLO. Partial NLO 
results on polarization can be found in [57, 58]. In this section we concentrate on one 
aspect of the problem, namely on the so-called anomalous contribution that, in the chiral 
limit, flips the helicity of the final-state lepton at NLO. 

Collinear photon emission from a massless fermion line can flip the helicity of the 
massless fermion contrary to naive expectation. This has been discussed in a variety of 
physical contexts. This is a "m e /m e " effect where the m e in the numerator is a spin flip 
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factor and the m e in the denominator arises from the collinear configuration. In the limit 
m e — > the helicity flip contribution survives whereas it is not seen in massless QED. 

We shall discuss this phenomenon in the context of the left-chiral [i — > e transition. At 
the Born term level an electron emerging from a weak (V — A) vertex is purely left-handed 
in the limit m e = 0. Naively, one would expect this to be true also at 0(a) or at any order 
in a because in massless QED photon emission from the electron is helicity conserving. 
Let us make this statement more precise by looking at the string of 7-matrices between the 
initial state [i spinor and final state e antispinor of the left-chiral /i — > e transition. Using 
(1 — 75) = (1 — 75) (1 — 75) /2, and the fact that in massless QED every photon emission 
brings in two 7-factors (one from the vertex and one from the fermion propagator) one 
finds 

tz e r("V (a - l5 ) u , = ^u e {t + 75 )r (? V (i - 75)%, (56) 

where stands for the (even) 7 matrix string brought in by the emission of n photons. 
We have commuted the left-chiral factor (1 — 75) to the left end where it then projects 
out the helicity state A e = —1/2 from the electron antispinor thus proving the above 
assertion. 

Let us take a closer look at the anomalous helicity flip contribution in leptonic fi — > e 
decays by considering the unnormalized density matrix element p ++ of the final state 
electron which is obtained by setting cos 9 = 1 in (38) (remember that G 5 vanishes for 
m e — > and G 6 = in the Standard Model). One has 

-— — = -f3x r ((Gi + G 3 ) + (G 2 + G 4 )Pcos# P ). (57) 

Contrary to naive expectations one finds non-vanishing right-handed (+) contribu- 
tions which survive the m e — > limit when one takes the m e — > limit of the NLO 
contributions to (57) [3]. In fact, one finds 

, dT , {+ \ = T^ r o(f(l " *) 2 (5 - 2x)l - [(1 - x)\l + 2x)]Pcose P ). (58) 
dxdcosUp 127T VL J L J ' 

The result is rather simple. In particular, it does not contain any logarithms or 
dilogarithms. The simplicity of the right-handed contribution becomes manifest in the 
equivalent particle description of /x-decay where, in the peaking approximation, /z-decay 
is described by the two-stage process n~ — > e~ followed by the branching process e~ — > 
e~ + 7 characterized by universal splitting functions D n f/hf(z) [56]. The symbols nf and 
hf stand for a helicity non-flip and helicity flip of the helicity of the electron. In the 
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splitting process z is the fractional energy of the emitted photon. The off-shell electron 
in the propagator is replaced by an equivalent on-shell electron in the intermediate state. 
Since the helicity flip contribution arises entirely from the collinear configuration it can 
be calculated in its entirety using the equivalent particle description. 

The helicity flip splitting function is given by D h f(z) = az/(2n), where z = k /E' = 
(£" — E)/E' = 1 — x/x', and where k is the energy of the emitted photon. E' and E 
denote the energies of the initial and final electron in the splitting process. The helicity flip 
splitting function has to be folded with the appropriate m e = Born term contribution. 
The lower limit of the folding integration is determined by the soft photon point where 
E' = E. The upper limit is determined by the maximal energy of the initial electron 
E' = m^/2. One obtains 



dr« a f 1 , , 1 dT BoTn ' ( -\x') ,„ x, 

= — ax (f ) 

dxdcosOp 2n Jx x' dx'dcosOp x' 

= ^-T f 1 dx'(x' - x) ((3 - 2x) + (1 - 2x')P cos 9 P ) 

Z7T Jx ^ ' 

= -^-YMl-xf(h-2x) -(l-xf (1 + 2x) P cos 9 P ), (59) 
12n ^ ' 

which exactly reproduces the result (58). 

Numerically, the flip spectrum function is rather small compared to the 0(a s ) no-flip 
spectrum function. However, when averaging over the spectrum the ratio of the 0(a s ) 
flip and no-flip contributions amounts to a non-negligible (—12%), due to cancellation 
effects in the 0(a s ) no-flip contribution. 



13 Summary and Concluding Remarks 

We have discussed NLO corrections to a multitude of polarization observables in different 
processes including nonzero mass effects. The results are available in compact analytical 
form. They are ready for use in physics simulation programs which are reliable even at 
corners of phase space where mass effects become important. Present and planned exper- 
iments (TEVATRON Run2, LHC, BELLE, BABAR, r-charm factories) will be sensitive 
to these NLO Standard Model effects. Standard Model NLO corrections to polarization 
observables are needed as background for possible new physics contributions. Last but 
not least the NLO results are needed for NLO sum rule analysis' involving polarization 
observables. All calculations are based on the same one-loop and tree-graph matrix 
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elements. Results on rates agree with previous calculations. Results on polarization 
observables agree with previous calculations where available. We have checked various 
limits and found agreement with previous mass zero calculations, the Goldstone boson 
equivalence theorem and the equivalent particle description of the anomalous helicity flip 
contribution. Because of the various checks and the fact that we have essentially used one 
set of matrix elements as input to the calculations we feel quite confident that our results 
are correct. 
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